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<N . Abstract 

The energy dependence of forward pion double charge exchange reactions on light 
^ ■ nuclei is studied for both the Ground State transition and the Double-Isobaric- Analog- 

l/"") . State transitions. A common characteristic of these double reactions is a resonance-like 

| peak around 50 MeV pion lab energy. This peak arises naturally in a two-step process 

in the conventional pion-nucleon system with proper handling of nuclear structure and 
pion distortion. A comparison among the results of different nuclear structure models 
demonstrates the effects of configuration mixing. The angular distribution is used to 
| fix the single particle wave function. 

^ ■ PACS: 25.80Gn, 24.10Ht,21.60Fw 

<-> '. 

1 Introduction 

> ■ 

The energy dependence of the forward angle cross section for pion-nucleus double charge 
exchange (DCX) reactions has been shown to contain important information on the fun- 
damental pion-nucleon charge exchange process and has been investigated extensively both 
experimentally [1, 2, 3] and theoretically [4, 5]. See Ref. [6] for a review of the field in 1993. 

The systematic data on the energy dependence of the DCX reaction consistently show a 
peak at low energy, which exists for both the double isobaric analog state (DIAS) and the 
ground state (GS) transitions. There have been controversial theoretical explanations of this 
peak. While some authors [4] have claimed that the peak is evidence for the existence of a 
dibaryon, it has been shown [5] that it can be understood based on a two-step sequential 
process which involves solely conventional degrees of freedom of pions and nucleons. In Ref. 
[5], a study on 42_48 Ca and 14 C showed that the appearance of this peak at low energy 
is naturally explained by the distortion of the pion scattering wave without invoking the 
existence of a dibaryon. 

The purpose of the present study is to investigate the question of what information about 
nuclear structure can be obtained from the energy dependence of the forward DCX reaction, 
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more specifically, what is the effect of configuration mixing on the energy dependence. In the 
present work we demonstrate the relevance of configuration mixing through a comparison of 
different nuclear models. In particular we use the shell model and models of extreme L-S 
and extreme j-j coupling. 

Early work on 16 indicated the importance of nuclear structure in the DCX reaction[7]. 
In that case the calculation was based on the plane wave approximation for the pion wave 
function. One understands now that the distortion of the pion wave (by an optical potential) 
plays a very important role in the DCX reaction. For a quantitative understanding of the 
effect of nuclear structure on DCX reactions a suitable optical model is indispensable. 

The present work uses the theoretical framework of the sequential process, as was done in 
[5]. For the pion wave distortion we use the optical model[8] that originated from Ref. [9, 10], 
as extended to the multi-j case. 

The article is organized as follows: The second section discusses the GS transition 
16 0(7r+,7r-) 16 Ne and 40 Ca(7r+, 7r^) 40 Ti . The third section studies the DIAS transition 
18 0(7r + , 7r~) 18 Ne . The fourth section contains conclusions and discussion. 

2 Ground State Transitions: 16 and 40 Ca 
2.1 Nuclear structure 

For the GS transitions 16 0(7r + , 7r _ ) 16 Ne and 40 Ca(7r + , 7r~) 40 Ti , the initial states can be as- 
sumed to be closed shell states so that one has to deal with configuration mixing only in 
the final states. The nuclear structure of the two final state nuclei has been investigated by 
several authors [11, 12, 13, 14]. Here we apply those results and compare them to the results 
derived from the models of extreme L-S and j-j coupling. 

The SU(4) £g> SU(3) model is a limiting case of L-S coupling, where SU(4) stands for 
the Wigner supermultiplet SU(4) symmetry [15] which assumes the degeneracy of all four 
spin-isospin states of the nucleon, and SU(3) stands for the Elliott model [16] which deals 
with the (orbital) rotational symmetry. A combination of these two symmetries, i.e. the 
SU(4) ® SU(3) model (denoted as SU(3) hereafter), works well in the p-shell and fairly well 
in the ds shell [17]. 

Another limiting case of nuclear structure is the single particle shell model (SPSM) which 
puts nucleons into j-orbits solely according to the single particle energy, neglecting the two- 
body couplings. This is an extreme case of j-j coupling. The shell model lies between these 
two extremes as far as the single particle occupancy in orbits is concerned [14]. In general 
the shell model is more realistic than the SU(3) and SPSM models. 

For the DCX reaction on 16 the initial and final states are the ground state of 16 
and 16 Ne, respectively. The simplest ground state of 16 is a closed shell. For a better 
description of the excited states of 16 0, a mixing of ds orbits was introduced [11]. However, 
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this ds mixing will contribute mainly to the excited states, whereas our interest is only in 
the ground state of 16 0. Therefore, for the study of the DCX reaction on 16 0, we restrict 
ourselves to the closed shell state. 

For 16 Ne, we use the weak coupling approximation [18], i.e. the ground state of 16 Ne is 
taken as consisting of two protons in the ds shell plus two neutron holes in the p shell. We 
make the further simplification that only the angular momentum zero pairs of the proton 
and neutron sector are considered. 

In Ref. [7] a truncated pi-ds-si space was used. Due to the large difference in the 
shapes of the energy dependence (see Section 3) of forward DCX amplitude among different 
configurations we pursue a description in the full p- and ds-shell. 

We begin with the SU(3) scheme where the ground state of the proton sector can be 
written as, 

|gs,p>su3=|[H]01;(40)0;0>, (1) 

in which the symbol at the right hand side, |[/]ST; (A/i)L; J >, contains the following quan- 
tum numbers: [f] and (A/z) are the SU(4) and SU(3) irreps, respectively; ST represent the 
spin and isospin of the two nucleon state, L the orbital angular momentum, and J the total 
angular momentum. The ground (hole) state of the neutron sector is given by 

|gs,n>su3=|[H]01;(02)0;0>. (2) 

The above SU(3) representation can be transformed into the j-j coupling formalism as follows, 



gS,P>SU3 = 




|gs,n>su3 = v / Tl(P|)~ 2 ' 0> +\/il(P|)" 2 ' 0> ' ( 3 ) 

where, for example, |(si) 2 ; > and Kpi) -2 ; > stand for a state with two particles in si 
orbit and a state with two holes in the pi orbit, respectively; and both states have angular 
momentum zero. 

The shell model describes a more "realistic" case, and normally contains a mixture in 
the SU(3) irreps. We take the occupancy of single particle orbits from previous studies[ll, 
12, 13, 14]. The probabilities of the hole state in the piand pi orbits are 0.83 and 0.17, 
respectively. Therefore, one can write the shell model hole state as 

|gs,ra > SM = 0.91|(pi)- 2 ;0 > +0.41 1 (p f )" 2 ; >, (4) 

where the phase is determined by requiring that the L = term be dominant. In a similar 
manner the proton sector shell model wave function is determined as 

| gs ,p > SM = 0.32|(si) 2 ;0 > +0.32| (da) 2 ; > +0.89|(d|) 2 ; > . (5) 



3 



Another limiting case of the nuclear structure is the single particle shell model, in which 
the two protons occupy the ds and the two neutron holes occupy the pi orbit. 

The amplitudes of proton pairs in ds-shell orbits and the amplitudes of neutron hole pairs 
in p-shell orbits are listed in Table 1, where the last three rows give the amplitudes for the 
SU(3) model, the shell model (SM) and the single particle shell model (SPSM), respectively. 



Table 1: Single particle distribution over orbits 
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0.517 
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0.517 


SM 
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0.319 
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0.910 


0.415 


0.224 
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SPSM 








1.0 


1.0 








1.0 






The ground state of 40 Ca is taken as a closed shell state, whereas the ground state of 
40 Ti is considered as consisting of two proton particles in the f i shell and two neutron holes 
in the ds shell. The proton sector of 40 Ti is restricted to contain only one orbit, fi, thus 
one only needs to deal with the configuration mixing in the neutron sector. The structure 
of the neutron sector is approximated by that of the nucleus 38 Ca. In the SU(3) scheme the 
ground state of the neutron sector is 

|gs,n>su3= |[11]01;(04)0;0>, (6) 

where the notation is the same as for 16 Ne. 

For the shell model we carry out a schematic calculation similar to that for the structure 
of 18 [18], where the surface delta function interaction is assumed. The ground state of 
39 Ca, J n = | + is the single hole state of (da) -1 , whereas the excited state J 71 " = | + can be 

taken as the single hole excited state (si) -1 - There being no excitation of | + at low energy, 
we restrict ourselves to the truncated hole space (ds — si ) _1 . By imposing the requirement 
on the equality of the calculated and experimental binding energies[18], the ground state of 
38 Ca can be calculated as, 

|gs,n> SM = 0.975|(d f )- 2 ;0 > +0.222| (si )~ 2 ; >, (7) 

where the proton closed shell is neglected. 

In the single particle shell model, the two neutron holes of 40 Ti are in the d| orbit. The 
amplitudes of neutron hole pairs for 40 Ti are also listed in Table 1. The proton sector is not 
given since there is only one orbit, fz. 
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Figure 1: Angular distribution of the reaction 16 0(7r + , 7r~) 16 Ne at a pion laboratory energy 
of 65 MeV. The solid line represents the result of shell model, the dashed line gives the result 
of the SU(3) model and the dashed-dotted line the result of single particle shell model. 

It can be seen from the results of both 16 Ne and 40 Ti that, as far as the amplitude of 
particle (hole) pairs in different orbits is concerned, the shell model provides an intermediate 
case between the SU(3) model and the single particle shell model, as was shown in a previous 
study [14]. 

2.2 Angular Distribution 

The calculation of the DCX reaction follows a procedure developed in a previous paper [8], 
which extends the work in [19] (AGGK) to a multi-j case. Here we give the outline of this 
method. The DCX operator (in coordinate space) is expressed as 

F 12 (fe, k') = [F (r 1: r2 ) + ^( ri) r 2 )ai ■ e x a 2 • e 2 ]T_(l)T_(2), (8) 




5 



where JF and T\ are related to non-spin-flip (NSF) and double-spin-flip (DSF) processes, 
respectively, and they are calculated numerically treating pion wave propagation, distortion 
and charge exchanges. The neglect of the single-spin-flip term is for simplicity and is thought 
to be a reasonable approximation [8]. In the multi-j case it is more convenient to express 
the DCX operator in pair form: 

F= Yl G J,{jij2},{j 3 j4}( b j 1 b j 2 ) J ' ( 5 i3 5 i4) J > ( 9 ) 

where {jij'2} means j± > ji- The two-body DCX matrix elements G j^j^^rsH} consist of 
NSF and DSF terms for which the formula can be found in Ref. [8]. 

In multipole form the DCX operator looks the same as in the AGGK paper, 

F = n J2 F L{nnUjsj,}(k, k'^b)^ ■ (b} 2 ~a u ) L / (2L + 1), (10) 

L,{jlj'2},{jzj4} 

however, the actual expression of the matrix elements Fl is more complex and can be derived 
from the pair-form matrix elements G j^j^ihji} ■ 

In calculating the angular distribution, the parameters of the optical potential are taken 
from Ref. [5] which presents the result of fitting the elastic scattering. Another important 
ingredient in the calculation is the single particle wave function, especially the wave function 
of the orbits which are involved in the DCX reaction. The single particle wave functions 
are calculated by using a Woods-Saxon type of potential, the parameters being so adjusted 
that both the separation energy of the proton and neutron and the rms proton radius match 
the experimental values. However, this is still not enough to determine the single particle 
wave function uniquely. We take the fitting of the angular distribution as an additional 
requirement, since it is sensitive to this wave function. 

The angular distribution is calculated for 16 0(7r + , 7r") 16 Ne at the pion laboratory energy 
of 65 MeV (Fig. 1) and for the DCX reaction 40 Ca(7r+, 7r-) 40 Ti at energies of 49 MeV and 
56 MeV(Fig. 2). Data are taken from Ref. [1]. For oxygen the agreement between theory 
and experiment is good except for the largest angle. For calcium the shape is approximately 
correct but the magnitude is too small, especially at 56 MeV. As we shall see these energies 
are on the sides of the peak so that slight changes in energy change the magnitude a great 
deal. 

2.3 Energy Dependence 

Now that all the parameters are fixed, the energy dependence of the forward amplitudes 
and cross sections can be calculated without any adjustable parameters. As a first step, we 
calculate the forward amplitude for all possible transitions between single particle orbits. A 
somewhat surprising feature is that there exist drastically different patterns in some of the 
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Figure 2: Angular distribution of 40 Ca(7r + , 7r~) 40 Ti at a pion laboratory energy of 49 MeV 
and 56 MeV. The solid line represents the result of shell model, the dashed line gives the 
result of the SU(3) model, the dashed-dotted line the result of single particle shell model 
and the dotted line the result of adding 5% of ds neutron hole state. 



energy dependence of forward amplitudes. For example, Fig. 3 gives the energy dependence 
of the forward amplitudes of DCX on 16 for the pair transitions of pi — > si and pi — > d.3, 

2 2 2 2 

where the symbol pi — > Si means that a pair of neutrons in the pi orbit is converted to 

2 2 2 

a pair of protons in the si orbit, etc. The dashed-dotted line, dashed line and solid line 
represent the real part, imaginary part and absolute value of amplitudes, respectively. One 
can easily notice that the transition pi — > si has a peak at low energy (around 50 MeV), 
whereas the transition pi — > d3 has a peak at resonance. The origin of this drastic difference 

2 2 

in pattern of the energy dependence of forward amplitude is not yet understood, but one 
can expect from this result that the energy dependence of the forward cross section will be 
sensitive to the configuration mixing in the nuclear structure. 

For DCX on 40 Ca the situation is somewhat different. In Fig. 4 four DCX amplitudes 
are given: si, da, ds and fV— >fa. While the first three transitions are involved in the DCX 

2 2 2 2 2 

on 40 Ca, the fourth, fa — >fa is given only for comparison. The first three transitions all have 
a very small value at resonance region. In contrast, the transition fa — >fa has a second peak 
at resonance. 
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The forward cross sections of DCX on 16 and 40 Ca are calculated with the configuration 
mixing ratio provided by the three nuclear models in Table 1. From table one can easily see 
the relevant transitions which are: pi, p3— > si, d3 and ds for 16 0, and si, d.3 and d 
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for 4U Ca, respectively The results of DCX forward scattering are given in Figures 5 and 6. 
The solid line, dashed line and dashed-dotted line represent results of the shell model, the 
SU(3) model and the single particle shell model, respectively. 
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Figure 3: Examples of the energy dependence of forward amplitudes for DCX on 16 0. Transi- 
tions are from the pi orbit to ds orbits. The dashed-dotted, dashed and solid lines represent 

2 

the real part, the imaginary part and the absolute value, respectively. Optical model pa- 
rameters are taken from fits to the elastic scattering. 



Figure 5 shows that with the shell model wave function (solid line) both peaks at low 
energy and at resonance are reproduced fairly well. While the SU(3) result (dashed line) is 
close to that of shell model, the single particle shell model gives a poor prediction. If one 
believes the optical model we use, then the energy dependence of the forward DCX reaction 
suggests that the SU(3) model is a good approximation for the wave function of 16 Ne. 

For DCX on 40 Ca, Fig. 6 demonstrates that the peak at low energy is reproduced by all 
three nuclear models but with different height. While the peak of shell model is lower than 



8 





Figure 4: Energy dependence of forward amplitudes for DCX on 40 Ca. Transitions are from 
ds orbits to the fi orbit. The dashed-dotted, dashed and solid lines represent the real part, 
the imaginary part and the absolute value, respectively. Optical model parameters are taken 
from fits to the elastic scattering. 

the experimental one, that of SU(3) is higher than the data. This suggests that our shell 
model calculation needs to be improved. A possibility that suggests itself is the addition of 
the configuration mixing of (ds) -2 . A mixing of this kind with 5% probability was tried, 
and the calculation gives a peak at low energy much closer to experiment. The angular 
distributions are also improved considerably. The improved results for angular distribution 
and forward cross section are shown with a dotted line in figures 2 and 5, respectively. This 
mixing of (d|)~ 2 seems reasonable, since its neglect in subsection 2.1 was only because of 
the absence the single hole state (d|) _1 in the spectrum of 39 Ca. However, due to pairing 
between neutrons, the ground state of 38 Ca may well contain a significant mixing of the 
component of |(ds) -2 ; >. However, the adding of the configuration ds for 40 Ti needs to be 
justified by independent nuclear structure studies. Here we take this as an example of the 
sensitivity of DCX amplitudes to configuration mixing. 

In both figures 5 and 6, although the height varies, the peaks appear at low energy for all 
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Figure 5: Energy dependence of forward DCX cross section on 16 0. For the meaning of the 
curves see the caption of Fig.l. 



nuclear models, which is consistent with previous claims that this peak is the consequence 
of the pion wave distortion. 

While the position and height of the peak around 50 MeV are reproduced reasonably well, 
there are two discrepancies concerning the second peak, i.e. the peak at around resonance. 
First the position of the calculated second peak for 16 O appears at around 180 MeV whereas 
the experimental one is at around 160 MeV. In a two-step sequential process (which is the 
framework of this paper) the second peak is produced by the resonant behavior of the p- 
wave SCX amplitude which corresponds to the delta-resonance although the explicit degree 
of freedom of the delta is not included. To shift the peak from 180 to 160 MeV one may 
need to explore the processes beyond the degrees of freedom of nucleons. 

The second discrepancy is that the calculation of 40 Ca shows almost no peak at resonance. 
Among the three transitions, si, ds, and ds^fi, which are involved in the DCX on 40 Ca, 
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Figure 6: Energy dependence of forward DCX cross section on 40 Ca. See the caption of Fig. 
2 for the meaning of the curves. 

the first two give almost no peak at the resonance (see Fig. 4), and the third one has a 
scarcely noticeable DCX amplitude in this region. This means that it is impossible to get a 
good fit in the resonance region by only adjusting the nuclear wave functions, in contrast to 
the case of 16 0, for which the magnitude of the DCX amplitudes at resonance is comparable 
to the low energy region. To reproduce the second peak on 40 Ca we may have to develop a 
better understanding of the energy dependence of DCX amplitude for transitions between 
orbits. 

Another possible remedy is to readjust the true absorption, since one of the consequence 
of true absorption is to reduce the amplitude at resonance. In this study all the optical 
model parameters, i.e. the ranges of six partial waves, the absorption and the s and p 
wave multipliers are taken without modification from the corresponding elastic 7r-nucleus 
scattering fit [20]. The readjustment of the true absorption requires a re-calculation of the 
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elastic 7r-nucleus scattering. 

Despite all these problems, we can say that the fair agreement between data and theoret- 
ical calculation, at least for low energies, is an indication of the validity of the optical model 
employed. 

3 Double Isobaric Analog State Transition 

The DCX reaction 18 0(7r + , 7r~) 18 Ne has been investigated experimentally by a number of 
authors [21, 22] with an emphasis on pion energies above 100 MeV, leaving the peak at low 
energy poorly explored. 

The structure of 18 and 18 Ne has been studied based on the analysis of spectra and 
transitions [12, 22], and it has been shown that the ground states of 18 and 18 Ne contain 
a significant mixing between two particle and four-particle-two-hole configurations. The 
ground state of 18 can be written as [12], 

|0+ >= 0.842|(nd|) 2 > +0.440|(nsi) 2 > -0.313|^ 00 >, (11) 

where the \^/ 00 > is a deformed state formed by exciting a pair of protons from the pi shell to 
the ds shell. The |\I/oo > wave function was given in Ref. [12] as a mixture of different four- 
particle-two-hole states. For the sake of simplicity, we consider only the angular momentum 
zero pairs while keeping the total wave function \^ 00 > normalized. Therefore, in this 
approximation, 

l*oo > = {^\( P ds finds) 2 > +A|( p d f ) 2 (nd | ) 2 >+^|(ps.) 2 (ns | ) 2 > 
+ I ^=|(pd f ) 2 (nd | ) 2 + (nd f ) 2 (pd3) 2 > 

+^|( P d|) 2 (nsi) 2 + (nd f ) 2 (^) 2 > 

+l|( P d|) 2 (^) 2 + (nda) 2 (^) 2 >}|G^r 2 >, (12) 

where the superscript 2 (-2) represents a two-particle (two-hole) state, and the angular 
momenta of all the pairs are taken to be zero. The ground state of 18 Ne is assumed to be 
the same as 18 except that neutrons are changed to protons and vice versa. 

We take the optical model parameters from 16 [20]. As was done in the case of the 
GS transition on 16 and 40 Ca we fit the angular distribution in order to determine the 
single particle wave function. In the case of the DIAS transition on 18 one has to pay more 
attention to the single particle wave functions for the following reasons. 

The structure of 18 is quite similar to that of 42 Ca: two neutrons outside of the doubly 
closed shell. In a rough approximation both nuclei can be described as a two-particle state. 
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Nuclear structure studies [12, 23] have shown that they both have a mixing between two 
particle states and four-particle-two-hole states, even the mixing rates are quite similar one 
to another. It has been found for 42 Ca [24] that the radius of the valence neutrons is much 
larger than what one would expect from a calculation with a Woods-Saxon type of potential. 
To be more specific, a Woods-Saxon potential calculation for 42 Ca gives an rms radius of 
the fi valence neutron of 3.95 fm, whereas the analysis in Ref. [24] of pion scattering data 
gives a value for this radius of 4.95 fm. This one femtometer difference may be understood 
as a consequence of the fact that the two neutrons form a kind of cluster, which increases 
the attraction between the two valence neutrons and diminishes the interaction between the 
nucleon pair and the doubly closed core. Due to the similarity in nuclear structure between 
18 and 42 Ca one may suspect that the same phenomenon occurs in 18 0, i.e. an rms radius 
of valence neutrons considerably larger than that derived from a Woods-Saxon potential. 
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Figure 7: Angular distribution of the DCX of 18 0(tt+, 7r~) 18 Ne at the pion laboratory energy 
of 164 MeV. The solid line presents the results with a WS1 type of wave function, the dotted 
line shows the WS type of wave function and the dashed line corresponds to the HOI type 
of wave function. See the text for the definition of types of wave functions 
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Figure 8: Energy dependence of forward DCX cross section on 18 0. The dashed-dotted line, 
the dashed line, the solid line, and the dotted line correspond to the HO0,HO, HOI and 
H02 type of wave functions, respectively. See text for the definition of the wave function 
type. 



With this assumption in mind, the angular distribution at pion energy 164 MeV is cal- 
culated with three types of single particle wave functions. The first type is from a standard 
calculation with a Woods-Saxon type of potential as described in Section 2. The rms radii 
of the Id orbit (rid), the 2s orbit (r ld ) and the total neutron sector (r n ) are 3.68 fm, 3.91 
fm and 2.86 fm, respectively. This wave function is denoted as the WS wave function, and 
the angular distribution calculated with this wave function is given in Fig. 7 by the dotted 
line. The second type is made by combining the 16 core of the WS wave function and a 
valence wave function from an expanded Woods-Saxon potential well. Its radii are given in 
Table 2. This type is denoted as WS1 and the resulting cross section is shown in Fig. 7 by 
the solid line. A third wave function is made in a similar way as the WS1 type except that 
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the valence wave function is from a harmonic oscillator potential well. Its radii are given in 
Table 2 under the label HOI and the corresponding angular distribution is given in Fig. 7 
by the dashed line. The radii of the WS1 and HOI wave functions are chosen to be very 
close in order to explore the dependence of angular distribution on details of wave function 
with their radii being the same. We note that in all cases only the shell model nuclear wave 
function (described by Eqs. (11) and (12)) was used. The experimental data are taken from 
Ref. [21]. 

The WS and WS1 types of wave functions differ in their and r 2s radii. A comparison 
between them shows that with increasing radius the angle of the first minimum decreases, 
and the agreement with data is improved considerably. This can be viewed as a strong 
support for the assumption of the enlargement of the radius of the valence neutron. Radii 
of rid — 4.66 fm and r2. s = 4.42 fm would appear to be reasonable estimates for the valence 
nucleons. 

A comparison of the cases of HOI and WS1 shows that these two wave functions, with 
almost equal radii for valence neutron, give nearly the same angular distribution despite the 
differences in the shape of their wave functions. Therefore, one may expect that in the study 
of the DCX reaction (at least at this energy) the nuclear radius plays a crucial role while 
the details of the wave function are of little importance. 

To further explore the effects of the enlargement of valence neutron radius we calculated 
the energy dependence of the forward DCX with the following four types of wave function: 
HOO, HO, HOI and H02, where the HOI type was explained above and the HOO, HO and 
H02 are made in the same way as of the type of HOI but with different radii. The cross 
sections resulting from wave functions of the types HOO, HO, HOI and H02 are represented 
in Fig. 7 by dashed-dotted, dashed, solid and dotted line, respectively. We list the radii for 
all types of wave function in Table 2. We used the harmonic oscillator type of wave function 
for the valence neutron for convenience, and it should not give results significantly different 
from those with an expanded Woods-Saxon potential. 

A very interesting feature of Fig. 8 which shows the forward angle energy dependence is 
that the height of the peak at 50 MeV increases with a decreasing radius of valence neutron 
while the height of the second peak at resonance is little changed. The following could be a 
possible explanation: The angular correlation is being held fixed (i.e. the shell model wave 
function) so, as the radius is increased the two nucleons are getting farther apart. The DCX 
cross section drops rapidly as that happens (like the inverse of the distance between them 
squared). Thus, at low energy where the absorption effects are less important, the cross 
section rises as the radius is decreased. On the other hand, as the radius is increased both 
nucleons are farther from the nucleus and the probability of being absorbed is decreased. 
This effect is very strong in the resonance region so that around the delta the two effects 
compete and tend to cancel. 

The forward scattering shown in Fig. 8 seems to favor a smaller radii which would 
contradict what was suggested by the angular distribution. However, the above treatment is 
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only schematic, a more sophisticated self-consistent calculation is needed to understand the 
enlargement of the neutron radius of 18 O. 

Table 2: Radii of the valence neutron and the total neutron sector for 18 0. See the text for 
the meaning of the wave function type. 



wave function type 


rid 


r-2s 


r n 


WS 


3.68 


3.91 


2.86 


WS1 


4.66 


4.42 


3.14 


HOO 


3.22 


3.22 


2.75 


HO 


3.70 


3.70 


2.87 


HOI 


4.65 


4.65 


3.13 


H02 


5.69 


5.69 


3.46 



The calculation reproduces only the gross feature of experimental data. As in the case 
of the GS transition, a peak appears at low energy (around 50 MeV), but the height of the 
peak is much lower than the experiment if one takes the angular-distribution favored radii 
(rid — 4.66 and r 2s = 4.42 fm). Another discrepancy is that the second peak is calculated 
to be around resonance, but the experimental one is around 120 MeV. 

In Ref. [21] an attempt was made to reproduce the peak at 120 MeV with a two- 
amplitude model, which mixes the DIAS amplitude of 18 and the GS transition amplitude 
of 16 0. Although it was only a schematic treatment, the reproduction of the second peak is 
an indication that a configuration mixing more comprehensive than what was done in this 
work might be successful. 

4 Conclusion and Discussion 

In summary, the present work has demonstrated that the conventional sequential process 
can reproduce fairly well the angular distribution and the energy dependence of forward 
DCX reactions on both the GS transition of 16 and 40 Ca and the DIAS transition on 18 0. 
In our calculation the peak at low energy arises naturally as the consequence of pion wave 
distortion without invoking explicit dibaryon degrees of freedom, 

It is often thought that the DCX reaction depends mainly on the reaction mechanism 
with detailed nuclear structure being unimportant. However, this work shows that nuclear 
structure is also important in describing the DCX reaction. Due to a drastic difference in 
the pattern of energy dependence of forward amplitudes among different transitions between 
single particle orbits, the DCX reaction is sensitive to the configuration mixing ratios in the 
nuclear structure. In this work three nuclear models were employed. While all three models 



16 



predict a peak at around 50 MeV, the height of the peak is quite different among them. The 
shell model gives fairly good predictions on the forward DCX reaction both for 16 and 40 Ca. 
The difference among the results of different nuclear models can give an indication of how 
large a role the nuclear structure plays in the DCX reactions. In other words, the energy 
dependence of the forward DCX reaction contains information not only on the fundamental 
charge exchange process but on the nuclear structure as well. 

In this work only the nucleon pairs with angular momentum zero were considered, i.e., 
only the elementary process (nn)^(7r + , 7r _ )(pp)^ was involved, where the superscript 
represents the angular momentum of the nucleon pair. We also made some calculations 
including the pairs with angular momenta 2 and 4, but no significant change in the DCX 
amplitude was found. However, the role of higher angular momentum pairs can not be 
excluded definitively. 

We looked at the consequences of enlarging the radii of orbitals in 18 O as might be 
expected from pairing of the valence neutrons. We found that such considerations can solve 
a significant fraction of the problem of the position of the first minimum but give worse 
agreement with the energy dependence of the forward cross section. 

The study in the present work can be extended to medium mass nuclei, as data show that 
the forward cross section of the DCX reactions 56 Fe(7r + , 7r~) 56 Ni and 93 Nb(7r + , 7r~) 93 Tc also 
has a peak around 50 MeV. These nuclei have a much larger shell model space which makes 
an exact shell model calculation almost impossible. However, some recent developments in 
nuclear structure theory can reduce the space substantially, which makes the study possible. 
One of the new models is the pseudo-SU(4) model [25] which deals with nuclei with nucleons 
in the pf-shell by using the technique of the SU(4) model while incorporating the large L-S 
coupling. We intend to extend the model to the mass region 28 < Z, iV < 50, and a study 
of the application of this model to DCX in this mass region is currently underway. 
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